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We report a measurement of the branching fractions of B — > D**£~ vt decays based on 417 fb _1 
of data collected at the T(4S) resonance with the BABAR detector at the PEP-II e + e~ storage rings. 
Events are selected by fully reconstructing one of the B mesons in a hadronic decay mode. A fit 
to the invariant mass differences m(D ( -* s> Tv) — m(D'*') is performed to extract the signal yields of 
the different D** states. We observe the B — > D**£~i?e decay modes corresponding to the four D** 
states predicted by Heavy Quark Symmetry with a significance greater than six standard deviations 
including systematic uncertainties. 
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PACS numbers: 13.20He,12.38.Qk,14.40Nd 



Scmileptonic B decays to orbitally-excited P-wave 
charm mesons (D**) are of interest for several reasons. 
Improved knowledge of the branching fractions for these 
decays is important to reduce the systematic uncer- 
tainty in the measurements of the Cabibbo-Kobayashi- 
Maskawa [l| matrix elements \V c b\ and \V u b\- For ex- 
ample, one of the leading sources of systematic uncer- 
tainty on \V c b\ measurements from B — ► D*£~p£ de- 
cays 0] is the limited knowledge of the background due 
to B ^ D**Tv t 0. 

The D** mesons contain one charm quark and one light 
quark with relative angular momentum L = 1. According 
to Heavy Quark Symmetry (HQS) [4j, they form one dou- 
blet of states with angular momentum j = s q + L = 3/2 
[Di (2420), 1)1(2460)] and another doublet with j = 1/2 
[£>o (2400), L>i (2430)], where s q is the light quark spin. 
Parity and angular momentum conservation constrain 
the decays allowed for each state. The D\ and D% states 
decay through a D-wave to D*tt and D^tt, respectively, 
and have small decay widths, while the Dq and D[ states 
decay through an S-wave to Dtt and D*ir and are very 
broad. 

B — ► D**£~Dg decays constitute a significant fraction 
of B semileptonic decays 0] and may help to explain the 
discrepancy between the inclusive B — > Xt~T>g rate and 
the sum of the measured exclusive decay rates 0, @, 0| ■ 
The measured decay properties for B — > D**t~vt can be 
compared with the predictions of the Heavy Quark Ef- 
fective Theory (HQET) Q. QCD sum rules [9] imply the 
strong dominance of B decays to the narrow D** states 
over those to the wide ones, while some experimental 
data show the opposite trend 13, Til ] . 

In this letter, we present the observation of B semilep- 
tonic decays into the four excited D mesons predicted 
by HQS and measure the B(B — > D**£~Pg) branching 
fractions. The analysis is based on data collected with 
the BABAR detector 12| at the PEP-II asymmetric-energy 



e + e~ storage rings at SLAC. The data consist of a total of 
417 fb _1 recorded at the T(4S) resonance, correspond- 
ing to approximately 460 million BB pairs. An addi- 
tional 40 fb — taken at a center-of-mass (CM) energy 40 
MeV below the T(4S) resonance, is used to study back- 
ground from e + e~ —►//(/ = u,d, s,c,t) continuum 
events. A detailed GEANT4-based Monte Carlo (MC) 
simulation [l3T ] of BB and continuum events is used to 
study the detector response, its acceptance, and to val- 
idate the analysis techniques. The simulation describes 
B D**£-v t _decays using the ISGW2 model Q, and 
non-resonant B — ► D^ir£~&e decays using the model of 
Goity and Roberts [lj|. 

We select semileptonic B — > D**£"vi decays with 
I = e, ji in events containing a fully reconstructed B me- 
son (-Btag), which allows us to constrain the kinematics, 



reduce the combinatorial background, and determine the 
charge and flavor of the signal B meson. D** mesons are 
reconstructed in the D^tt^ decay modes and the differ- 
ent D** states are identified by a fit to the invariant mass 
differences m(D^n) — m(D^). 

We first reconstruct the semileptonic B decay, select- 
ing a lepton with momentum p| in the CM frame larger 
than 0.6 GeV/c. We search for pairs of oppositely- 
charged tracks that form a vertex and remove those with 
an invariant mass consistent with a photon conversion 
or a 7T° Dalitz decay. Candidate D° mesons that have 
the correct charge correlation with the lepton are recon- 
structed in the K~ir + , K~ir + ir°, K~Tr + ir + ir~ , KgTr + Tr~ , 
K°ir+ir~ir°, K°ir°, K+R-, tt+tt", and K%K% channels, 
and D + mesons in the K~ir + Tr + , K~ir + ir + ir° , Kgir + , 
J^tt+tt , K+K-JT+, K° S K+, and i^Tr+Tr+vr- channels. 
In events with multiple D£~ combinations, the candi- 
date with the best D-£ vertex fit is selected. Candidate 
D* mesons are reconstructed by combining a D candi- 
date with a pion or a photon in the D* + — > D°ir + , 
D*+ -> D+ir°, D*° -> D°ir°, and D*° -> D°j channels. 
In events with multiple D*£~ combinations, we choose 
the candidate with the smallest x 2 based on the devia- 
tions from the nominal values of the D invariant mass 
and the invariant mass difference between the D* and 
the D, using the resolution measured in each mode. 

We reconstruct B tag decays [l6| in charmed hadronic 
modes B — > DY, where Y represents a collection of 
hadrons, composed of niit^ +n2K^ +713X5+714^° , where 
ni+n 2 = 1, 3, 5, n 3 < 2, and n 4 < 2. Using D°(D+) and 
D*°(D* + ) as seeds for B~(B°) decays, we reconstruct 
about 1000 different decay chains. 

The kinematic consistency of a -Btag candidate with 
a B meson decay is evaluated using two variables: the 
beam-energy substituted mass ities = \Z S /^ — bsl 2 , 
and the energy difference AE = E* B — ^/s/2. Here ^/s is 
the total CM energy, and p* B and E* B denote the momen- 
tum and energy of the -Btag candidate in the CM frame. 
For correctly identified B tag decays, the thes distribu- 
tion peaks at the B meson mass, while AE is consistent 
with zero. We select B tag candidates in the signal region 
defined as 5.27 GeV/c 2 < m ES < 5.29 GeV/c 2 , excluding 
those with daughter particles in common with the charm 
meson or the lepton from the semileptonic B decay. In 
the case of multiple B tag candidates in an event, we select 
the one with the smallest |A-B| value. The B tag and the 
D(*>£ candidates are required to have the correct charge- 
flavor correlation. We account for mixing effects in the 
B° sample as described in Ref. [IT] - Cross-feed effects, 
i.e., B t ~ g (B° ag ) candidates erroneously reconstructed as 
a neutral (charged) B, are subtracted using estimates 
from the simulation. 

We reconstruct B~ -> B ( *^i\-£rv t and B° -> 
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Mode 



Selection Criteria 



B~ 
B~ 
B° 
B° 



D* 



7T t 
TY~t 



D ^ 



-0.25 < m 
-0.25 < m 
-0.2 < m: 
-0.15 < m 



i — 

miss 
2 

miss 

2 

miss 
2 

miss 



< 0.25 GeV7c 4 

< 0.8 GeV 2 /c 4 
< 0.35 GeV 2 /c 4 

< 0.85 GeV 2 /c 4 



_D(*)°7t+£ vg decays starting from the corresponding 



B 



lag 



D(*H~ combinations. We select events with 
only one additional reconstructed charged track, cor- 
rectly matched to the D" flavor, that has not been used 
for the reconstruction of the -Btag, the signal D^*\ or the 
lepton. D(D*) candidates are selected within 2er (1.5- 
2.5cr, depending on the D* decay mode) of the D mass 
(D* — D mass difference), where the resolution a is typi- 
cally around 8 (1-7) MeV/c 2 . For the B° -► V+ i~ v t 
decay, we additionally require the invariant mass differ- 
ence m(D°n + ) - m(D°) to be greater than 0.18 GeV/c 2 
to veto B° — > D* + £~D e events. 

Semileptonic B — > D**£~vt decays are identi- 
fied by the missing mass squared in the event, 
< iss _ = [p{T{AS))-p{B tag )-p{D^^)-p{£)\\ de- 
fined in terms of the particle four-momenta. For correctly 
reconstructed signal events, the only missing particle is 
the neutrino, and rn^^ peaks at zero. Other B semilep- 
tonic decays, where one particle is not reconstructed 
(feed-down) or is erroneously added to the charm candi- 
date (feed- up), exhibit higher or lower values in 0- 
In feed-down cases where both a D and a D* candidate 
have been reconstructed, we keep only the latter candi- 
date. 

The m^iss selection criteria are listed in Table HI The 
m^jgg region between 0.2 and 1 GeV 2 /c 4 for B — > 
Dni^Dg events is dominated by feed-down from B — > 
D**(-^ D*ir)£~Di semileptonic decays where the soft 
pion from the D* decay is not reconstructed. In order 
to retain these events we apply an asymmetric cut on 
m^jgg for these modes. 

The signal yields for the B — > D**£~Dg decays are 
extracted through a simultaneous unbinned maximum 
likelihood fit to the four m(D ( -*K) - m(£>W) distribu- 
tions. With the current statistics, validation studies on 
MC samples show that our sensitivity to non-resonant 
B — ► D^*'TTi~9i decays is limited. Including hypothe- 
ses for these components results in a fitted contribution 
that is consistent with zero. Thus we assume that these 
non-resonant contributions are negligible. The probabil- 
ity that B — > D**(-^ D*tt)£~vi decays are reconstructed 
as B — > D**{-^ Dir)£~&£ is determined with the MC sim- 
ulation to be 26%(59%) for the B~(B°) sample and held 
fixed in the fit. 

The Probability Density Functions (PDFs) for the 
D** signal components are determined using MC B — > 
D**£~D( signal events. A convolution of a Breit-Wigner 




nD*tv 

I I background 



a) 



b) 



0) 



•i Hi 



d) 



0.8 ,„ 1 „, 1.2 

M(D ( , jc)-M(D () ) [GeV/c 2 ] 



FIG. 1: (Color online) Fit to the mfZ/'V) - m{D M ) dis- 
tribution for a) B~ -> D* + ir~t~vt, b) B~ -> D+tvTuh, c) 
B° -> D*°ir + rv e , and d) B° -> D°n + £~ v t : the data (points 
with error bars) are compared to the results of the overall fit 
(sum of the solid distributions). The PDFs for the different 
fit components are stacked and shown in different colors. 



function with a Gaussian, whose resolution is determined 
from the simulation, is used to model the D** resonances. 
The D** masses and widths are fixed to measured val- 
ues We rely on the MC prediction for the shape 
of the combinatorial and continuum background. A non- 
parametric KEYS function [18| is used to model this com- 
ponent for the D*n£~p£ sample, while for the Dit£~vi 
sample we use the convolution of an exponential with 
a Gaussian to model the tail from virtual D* mesons. 
The combinatorial and continuum background yields are 
estimated from data. We fit the hadronic i?tag m>ES dis- 
tributions for B — > D**£~V£ events as described in 0], 
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and we obtain the number of background events from 
the integral of the background function in the hues sig- 
nal region. 

Table [TT] summarizes the results from two fits: one in 
which we fit the charged and neutral B samples sepa- 
rately, and one in which we impose the isospin constraints 
B{B- -> D**e-v e )/B(B° -> D**l-vt) = t b -/t Bo . The 
latter fit yields a significance greater than 6 standard 
deviations for all four D** states including systematic 
uncertainties. The results of this fit are shown in Fig. 1. 

The D2 contributes to both the Dir and the D*ir 
samples. In the nominal fit we fix the ratio B(D\ — > 
Dtt)/B{DI -» D*tt) to 2.2 0. When we allow this ratio 
to float we obtain 1.9 ± 0.6. 

To reduce systematic uncertainties we measure the ra- 
tios of the B(B — > D**£~Pe) branching fractions to the 
inclusive B° and B~ semileptonic branching fractions. A 
sample of B — > X£~vi events is selected by identifying 
a charged lepton with p* t > 0.6 GeV/c and the correct 
charge correlation with the -Btag candidate. In the case of 
multiple Bt a g candidates in an event, we select the one re- 
constructed in the decay channel with the highest purity, 
defined as the fraction of signal events in the tues signal 
region. Background components that peak in the m B s 
signal region include cascade B meson decays (i.e., the 
lepton does not come directly from the B) and hadronic 
decays, and are subtracted using the corresponding MC 
predictions. 

The total yield for the inclusive B — > X£~p£ decays is 
obtained from a maximum likelihood fit to the rriES dis- 
tribution of the Stag candidates, as described in The 
fit yields 198,897 ± 1,578 events for the B~ -> Xl~v t 
sample and 120,168 ± 1,036 events for the B° -> Xl~v t 
sample. 

The ratios B{B -► D**£-D e )/B(B -> Xfvt) = 
(A^sig/csig) • (e s \/N s \) are obtained by correcting the sig- 
nal yields for the reconstruction efficiencies (estimated 
from BB MC events). Here, 7V s i g is the number of 
B — > D**i~Di signal events, reported in Table [TT] to- 
gether with the corresponding reconstruction efficien- 
cies e s i g , N s \ is the B — > X£~T>i signal yield, and e s i is 
the corresponding reconstruction efficiency including the 
Stag reconstruction, equal to 0.39% and 0.25% for the 
B~ — > X£~p£ and B° — > X£~vt decays, respectively. 
The absolute branching fractions B(B — > D**£~Dg) are 
then determined using the semileptonic branching frac- 
tion B{B -> X£~vi) = (10.78 ± 0.18)% and the ratio of 
the B° and the B~ lifetimes t b -/t b „ = 1.071±0.009 @]. 

Numerous sources of systematic uncertainties have 
been investigated. The largest uncertainty is due to the 
determination of the B — ► D**£~Dg signal yields (result- 
ing in 5.5-17.0% relative systematic uncertainty depend- 
ing on the D** state). This uncertainty is estimated using 
ensembles of fits to the data in which the input parame- 
ters are varied within the known uncertainties in the PDF 
parameterization (0.2-8.7%), the shape and yield of the 



combinatorial and continuum background (0.2-10.4%), 
the modeling of the broad D** states (4.5-13.8%), and 
the D* feed-down rate (0.5-4.0%). We check that the 
combinatorial and continuum background shape is well 
reproduced by the simulation by verifying that the MC 
samples of right-sign and wrong-sign D^n combinations 
have similar shapes, and that the wrong-sign distribution 
in the data agrees well with that in the simulation. Wc 
observe an excess of events in the low invariant mass dif- 
ference region in the four samples that is not accounted 
for by the background PDF. We study B — > D^mr£~ Dg 
(n > 1) decays, not included in our standard MC simula- 
tion, as a possible source of this excess. We use different 
MC models for these decays, and find that they do not 
account for all the observed excess. We evaluate a corre- 
sponding systematic uncertainty (0.1-3.2%), included in 
the yield uncertainty above. The uncertainties due to the 
detector simulation are determined by varying, within 
bounds given by data control samples, the charged track 
reconstruction efficiency (1.3-2.0%), the photon recon- 
struction efficiency (0.2-4.8%), the lepton identification 
efficiency (1.2-1.6%), and the reconstruction efficiency for 
low momentum charged (1.2%) and neutral pions (1.3%). 
We use an HQET model [8j to test the model dependence 
of the B -> D**£~v t simulation (0.8-2.5%). We include 
the uncertainty on the branching fractions of the recon- 
structed D and D* modes (3.0-4.5%), and on the absolute 
branching fraction B(B — > X£~vg) used for the normal- 
ization (1.9%). We also include a systematic uncertainty 
due to differences in the efficiency of the Stag selection in 
the exclusive selection of B — ► D**£~p£ decays and the 
inclusive B — > X£~Pi reconstruction (4.0-5.6%). 

In conclusion, we report the simultaneous observation 
of B — > D**£"Dg decays into the four D** states pre- 
dicted by HQS. The measured branching fractions are re- 
ported in Table HT1 We find results consistent with Ref. Q 
for the sum of the different D** branching fractions. The 
rate for the D** narrow states is in good agreement with 
recent measurements [li| ; the one for the broad states is 
in agreement with DELPHI [ll[ but does not agree with 
the D[ limit of Belle [10( ■ The rate for the broad states is 
found to be large. If these broad states are indeed due to 
B — ► D[£~Di and B — * Dq£~Di decays, this is in conflict 
with the expectations from QCD sum rules. 
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viduals have received support from the Marie Curie EIF 
(European Union) and the A. P. Sloan Foundation. 
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TABLE II: Results from the fits to data: the B — > D**t~vt signal yield, the corresponding reconstruction efficiency, the product 
of branching fractions, where the first error is statistical and the second systematic. For the B — > D^l vg decay, we report 
yields and product of branching fractions for the D| — > Dn decay mode. For the isospin-constrained results (last two columns), 
the B~ branching fraction products are reported. The statistical significances, S'stat, are obtained by computing the difference 
in the log likelihood between the nominal fit and the fit in which we fix the different signal components to 0. The significances 

including the systematic uncertainty, S to t, are obtained by rescaling the statistical significances by CstatA/oftat + °"?ysf 



Decay Mode 


Yield 


esig (xl(T 4 ) B (B-> 


D**Tvi ) x B(D" -> D^ ) n ± ) 


% S-tot {Sstat ) 


B % Stot (iSstat) 


B~ 


Dirvi 


165 ± 18 


1.24 




0.29 ±0.03 ±0.03 


9.9 (12.7) 


0.29 ±0.03 ±0.03 10.7 (15.2) 


B~ 


— Dfrvt 


97 ± 16 


1.44 




0.15 ±0.02 ±0.01 


6.3 (7.3) 


0.12 ±0.02 ±0.01 6.0 (7.4) 


B~ 


-> D^rvt 


142 ± 21 


1.13 




0.27 ±0.04 ±0.05 


5.4 (8.0) 


0.30 ±0.03 ± 0.04 6.4 (10.0) 


B~ 


-> Dfrut 


137 ± 26 


1.15 




0.26 ± 0.05 ± 0.04 


4.5 (5.8) 


0.32 ±0.04 ± 0.04 6.1 (8.3) 


B° 


-> D+Tvt 


88 ± 14 


0.70 




0.27 ±0.04 ±0.03 


7.0 (8.4) 




B° 


-» D*+i-9 t 


29 ±13 


0.91 


0.07 ± 0.03 ± 0.01 (< 0.11 @90% CL) 


2.3 (2.5) 




B° 


-> D+rv t 


86 ± 18 


0.60 




0.31 ±0.07 ±0.05 


4.6 (5.8) 




B° 




142 ± 26 


0.70 




0.44 ±0.08 ±0.06 


4.7 (6.0) 
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